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Abstract 

We discuss some of the challenges that future nuclear modeling may face in 
order to improve the description of the nuclear structure. One challenge is related 
to the need for A-body nuclear interactions justified by various contemporary nu- 
clear physics studies. Another challenge is related to the discrepancy in the NNN 
contact interaction parameters for 3 He and 3 H that suggests the need for accurate 
proton and neutron masses in the future precision calculations. 

1 Introduction 

The high precision, QCD derived, nucleon interaction that describes the NN-scatering 
phase shifts, the deuteron, and the light s- and p-shell nuclei points to the necessity of 
NNN-interaction terms |T][2j. Thus the conventional two-body interaction paradigm 
is challenged and the need of 3-body and possibly A-body interaction define a new 
research frontier. The structure of the three-body terms has been studied previously 
using the meson exchange theory (3). However, with the advance of the Chiral Per- 
turbation Theory (ChPT) [4 5| the structure of the three-body terms is better justified 
using QCD. While studying the parameters of the 3-body contact terms (2j one faces 
a discrepancy in the NNN contact interaction parameters needed to fit 3 He and 3 H that 
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could be viewed as an argument towards implementing the accurate proton and neutron 
masses in future precision calculations. 

Higher many-body interaction terms (e.g. NNNN-interaction terms) are also part 
of the interaction as derived from QCD via ChPT |6|. The Okubo-Lee-Suzuki (OLS) 
effective interaction method, employed in solving the nuclear many-body theory, also 
introduces interaction terms beyond the common 2-body interaction |7][8). All this 
seems to be pointing to the need of A-body interactions for the description of the nu- 
clear structure. It also raises the question about the importance of the A-body inter- 
actions in very heavy nuclei. Fortunately, there is an exactly solvable A-body model - 
the extended pairing model - that is applicable as an A-body interaction to very heavy 
nuclei; therefore, it can help to address this question ]9]-|TT) . 

In the next section we briefly discuss the microscopic nuclear physics hamiltonian; 
the types of the high-precision NN-interaction potentials and their failure to properly 
account for the structure of the nuclei with more than two nucleons. In Sec. [3] we 
discuss the discrepancy in the values of the cjj and <?£ NNN-intercation parameters |2| 
and try to argue that the application of high-presision nucleon potentials needs more 
appropriate nucleon masses for high-precision description of the A=3 systems. In Sec. 
[4]we further extend our argument for A-body nuclear interactions by using the modern 
OLS effective interaction in finite model space method. In Sec. |5]we briefly discuss the 
well-know 2-body pairing interaction and its exact solution as a prelude to the A-body 
Extended Pairing Interaction (EPI); then we discuss the results of apply the EPI to few 
long isotope chains like Sn, Yb, and Pb nuclei. Last section is our conclusion about the 
needs of the future nuclear structure modeling. 

2 Modeling the Nuclear Interactions 

Unlike the electromagnetic and the gravitational intersection, the mathematical form 
of the nuclear interaction has been very elusive. It is now clearly understood that 
this is due to the fact that the nuclear interaction arises nontrivially from the quark 
structure of the nucleons and thus related to the theory of the QCD. However, the 
absence of a closed form interaction has not hindered researchers from modeling the 
structure of nuclei. The field has advanced significantly, based on general quantum 
mechanical principals and techniques. In particular, the microscopic approach has been 
very successful especially with the advance of computational techniques and computer 
power that have allowed for the construction of effective high-precision meson and/or 
QCD derived NN-potentials. The free parameters of the high-presision NN-potentials 
are usually fixed by the experimental two-nuceon scattering data and describe the 2- 
body system extremely well. Unfortunately, these potentials produce unsatisfactory 
description of the 3- and 4-body systems. 

2.1 The Nuclear Shell-Model Hamiltonian 

A nuclear many-body system near equilibrium can be viewed as subject to a mean field 

Harmonic Oscillator (HO) potential: #0 = 5^ + \k 2 x 2 . It is well know that one can un- 
derstand the magic numbers and the shell structure of nuclei within the 3-dimensional 
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HO approximation plus a spin-orbit potential [ 12 1. Using the HO single-particle states 
one can write a general Hamiltonian with one- and two-body terms: 

H = ^eiafai + - £ V^/a+a+a^a/. (1) 

i i,j,k,l 

Here, a, and aj are fermion annihilation and creation operators, £; single-particle ener- 
gies, and Vuki = (ij\V\kl) two-body interaction matrix elements and the index i labels 
the single particle levels. Despite the significant symmetry relations, e.g. Ej m = £,„,/ 
due to rotational symmetry and Vum — Vkijj — —Vnu — —Vijik due to the fermion 
exchange properties and the hermition requirement on the energy operator, the number 
of independent parameters is often more than a dozen - usually it is of order of few 
hundred for the valence NN interactions alone. The independent parameters of the in- 
teraction ([T| are often fitted to experimental data by starting with some initial values 
that come from a relevant theory or model. 



2.2 Problems with the High-Precision NN-Potentials 

Many of the high-precision NN-potentials commonly used to build the microscopic 
interactions for multi-nucleon systems have very complicated but methodically devel- 
oped structure in terms of spin, iso-spin, and angular momentum components although 
sometimes there is a very complicated radial dependence. For example, the Argonne 
V18 potential has 18 different terms |T3| . Other potentials use non-local terms e. g. 
CD-Bonn p4[ and Nijmegen | |T5| . However, when applied to A>2 systems all of these 
potentials have a serious difficulties that were eventually overcome by using three-body 
interactions (4][T6j[T7) . 

By the end of the twentieth century it become clear that a two-body interaction by 
itself is inadequate even for the description of the lightest nuclei 2 < A < 5. Com- 
parative studies of various potentials, such as AVI 8, Nijmegen, CD-Bonn, and N 3 LO, 
with or without three body terms have demonstrated the inadequacy of the pure two- 
body interactions and the need for three-body interaction terms p7][T8) . For example, 
all these interactions (AVI 8, Nijmegen, CD-Bonn, and N 3 LO) describe very well the 
deuteron properties such as binding energy, radius, and quadruple moment but fail by 
more than 0.5 MeV to reproduce the binding energy of triton [ 18] and underbind 4 He 
by more than 4 MeV fT7) . 

Although the meson-exchange approach was successful, it was clear that this phe- 
nomenological models should be derived from the underling QCD. Thus the ChPT ap- 
proach became a prominent technique that produced the high-precision NN-potential 
N 3 LO and then guided the researchers into the structure of the NNN- and NNNN- 
interactions (6][T8j[T9] • 



3 



3 Light Nuclei and the Parameters of the NNN-body 
Interaction 



The use of the ChPT in the derivation of the nucleon interactions from QCD assisted 
in the determination of the mathematical form of various interaction terms along with 
the relevant parameters. Unfortunately, parameters related to contact terms in the in- 
teraction could not be determined. Thus the cd and the ce strengths of the two-nucleon 
contact interaction with one-pion exchange to a third nucleon and the three-nucleon 
contact interaction are identified freedoms at the present time in the effective ChPT 
interaction. As such they need to be fixed by comparison with experiment. 

3.1 Binding Energy of 3 H, 3 He, and 4 He 

In order to determine the cp and ce parameters of the interaction one searches for 
the parameter values that reproduce the binding energy of 3 H and 3 He within 0.5 keV 
of the experimental values [2|. As seen from Figure [T] there are two cd — ce curves 
that unfortunately do not intersect. In order to further narrow down the range of cd 
values one considers the averaged cp — ce curve and evaluates the binding energy of 
the 4 He system which results in two possible physical regions denoted by A and B on 
inset (a) of Figure [T] Finally, the charge radius of 4 He points to the region A as the 
reasonable range of values for the cp parameter while the ce parameter is determined 
by the averaged co — ce curve. 
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Figure 1: Relations between co and ce for which the binding energy of 3 H (8.482 
MeV) and 3 He (7.718 MeV) are reproduced, (a) 4 He ground-state energy along the 
averaged co — ce curve. The experimental 4 He binding energy (28.296 MeV) is re- 
produced to within 0.5 MeV over the entire range depicted, (b) 4 He charge radius r c 
along the averaged co — ce curve. Dotted lines represent the r c uncertainty due to the 
uncertainties in the proton charge radius. 
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Conceptually, there are three important concerns: First, the ChPT NN-potential 
was one order higher than the NNN-potential and no NNNN-potential was included. 
That is, the high-precision NN-potential was N 3 LO (next-to-next-to-next-to-leading 
order) fllflw hile the ChPT NNN-potentail was at the N 2 LO order (l9) and the NNNN- 
potential Q was not yet readily available. The second concern is that the range of the 
3-body interaction parameter cu is determined by the properties of the 4-body system 
4 He; this, however, was resolved by a later study that used the j3~ decay of 3 H into 
3 He and confirmed the physically relevant region A for the parameter cd Q20J. The 
third concern is related to the fact that these are high-presision studies and at this level 
of accuracy the difference between the proton and nucleon mass could be important for 
the A=3 systems. 

3.2 3 H and 3 He Systems with Modified Nucleon Mass 

It is clear from Figure [T] that the co — ce curves for 3 H and 3 He do not intersect in 
the physically relevant region (—1 < co < 1). This could be attributed to the absence 
of the T=3/2 channel in these first calculations. The slight difference in the cp value 
as suggested by 4 He binding energy and its charge radius could be attributed to the 
inconsistency of the different interaction terms, i.e. NN-terms are at N 3 LO level while 
the NNN-terms are at N 2 LO level and the NNNN-terms are not present at all. Another 
source of these discrepancies could be the conventional use of equal masses for protons 
and neutrons. We will discuss this option in more detail below. 

As seen from Figure [5] the binding energy deviation range is 15 < SK < 25 keV 
within the physically relevant region (—1 < cq < 1). This is within the accuracy of the 
kinetic energy K as evaluated for equal mass nucleons m = m„ = m p = (m„ +m p )/2. 
Since the averaged relative kinetic energy for the three-nucleon system is about K « 37 
MeV and the relative nucleon mass deviation Sm/m is « 0.7 x 10~ 3 with 8m = (ra„ — 
m p )/2, we have: 

K= M v 2 ^ SK = K— «26keV (2) 
2 m 

This shows that there is not a single ce value that will result in perfect description of 
the 3 H and 3 He systems. One could hope that including the T=3/2 channel would im- 
prove the situation. Alternatively, with this level of precision, we are led to investigate 
corrections to the conventional m„ = m p approximation. One can test the sensitivity to 
the conventionally used nucleon mass by changing it to a more appropriate value pT) . 

m = i (Zm p + (A - Z)m„). (3) 

If one repeats the calculations related to Figure [2] but by employing the nucleon 
mass value as suggested by ([3]), one obtains an interesting result shown in Figure [3] 
From Figure [3] is clear that there are unique ce> and eg values where both binding 
energies can be reproduced exactly. The cd is in agreement with the cd value estimated 
from the charge radius of 4 He (see Figure[TJb)). Since, 4 He has same number of protons 
and neutrons, there is no nucleon mass adjustment for this system, so results shown for 
4 He on Figure [l] are still valid. Perhaps by incorporating NNNN-interaction |6| the 
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Figure 2: Deviation of the binding energy of the three-nucleon systems as computed 
along the averaged cd — ce curve. 

binding energy for the 4 He would agree better with the cp value suggested by the 
three-nucleon system as calculated with a modified nucleon mass and the j3 ~ decay of 
3 H (20). 

4 Beyond the 2-body Interaction - Effective Interactions 
in a Finite Model Space 

In the previous section we discussed results obtained by using QCD derived interac- 
tions and the role of the NNN-interaction in the description of the light nuclei. Clearly 
3- and 4-body interaction terms are predictions of the ChPT. Thus A-body interactions 
can be viewed as real physical interactions within the ChPT approach to nuclei. How- 
ever, there is another way to arrive at A-body interactions which are phenomenological 
effective interactions since they are related to our inability to handle interacting systems 
in infinite Hilbert spaces (22). Since the quality of a model is judged by its ability to 
reproduce the experimental data, as far as computational models are concerned, an A- 
body interaction which gives results that agree well with the data is physically relevant 
as well. 

In practice, we are computationally limited to a finite subspace of the infinite 
Hilbert space of the full quantum many-body problem. The subspace that we can 
access is defined by finite set of convenient many-body basis states. For a suitable 
choice of basis we hope to have good overlaps with low-lying physical states of the 
system under study. If we imagine the exact solutions are available for analysis and 
apply a unitary transformation to those eigenstates, we can produce a transformed set 
of solutions maximally overlaping with our chosen basis space. 

For example, one may be interested in the lowest two energy states of a system, as 
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Figure 3: Intersecting cp — ce curves where the binding energy of 3 H (8.482 Me V) and 
3 He (7.718 MeV) are reproduced when using modified nucleon mass as suggested by 

shown in Figure |4] left, but would like to have some unitarily transformed version of 
these states that have maximal overlap with the two basis states that define the plane 
of the page (Figure |4] right). By finding the relevant unitary transformation U, one 
can define an effective Hamiltonian that would have the lowest two states as desired. 
Then this effective Hamiltonian could be used in the calculations of more complicated 
multi-particle systems, i.e. one would find the unitarily transformed Hamiltonian that 
describes very well the low-energy states of a 2-body system in a mean field but within 
a Fock space that would be used later for an A-body system. Unfortunately, this trans- 
formation will turn any one- and two-body potential into a many-body effective inter- 
action: 

i A o A i A 

This way the two-body Hamiltonian ([TJ becomes an A-body Hamiltonian: 

A l 

H = ^e/fl+a; + £ —2 £ Viu---,i k ,h,---,h a h ■ ■■ a X a h ■ ' ' a k ■ ( 4 ) 
i k=2 y K -> 

ll,...,Jfc, 

jli—Jk 

For A>4, it seems impractical at present to obtain the structure of the A-body 
interactions as derived from ChPT as it was previously done for the NNN- and the 
NNNN-interaction terms. Before embarking on the extensive undertaking required for 
including higher-body effective interactions, it would be very helpful to investigate 
a simple exactly solvable A-body interaction model that has few parameters and is 
applicable to real A-body systems. 
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Figure 4: Geometrical interpretation of the Okubo-Lee-Suzuki transformation method 
for construction of effective Hamiltonian operators. 

5 The Extended Pairing Model 

In order to study the relevance of the A-body interactions one should use the general 
form of the interaction and to try to determine some of the A-body interaction strengths 
since it seems impractical at present to be able to obtain the structure of the A-body 
interactions from ChPT for A>4. Therefore, as we reasoned earlier one needs simple 
exactly solvable A-body interaction with few parameters that can be adjusted to the 
experimental data. Fortunately, there is such an interaction - the Extended Pairing 
Interaction (EPI) |9|. The discovery of this exactly solvable model was a result of 
research into the solution of the two-body proton-neutron pairing which turned to be 
exactly solvable as well (23]. However, for our purpose the justification, of the A- 
body EPI Hamiltonian, is the need for simplicity: thus, one can set all the unknown 

interaction strengths Vq I( . j, j k in (HI to be equal to single strength G and to consider 

only pairs of fermion particles bf — atat: 



i i.j k=2\ K -> ii.-hk 

Here n,- counts the number of pairs on the 2-th level; thus, the value 2 in front of the 
single particle energy £,-. If one considers a system of only one pair of particles then the 
k > 1 terms in |5]) disappear since their matrix elements are zero in the one-pair basis. 
Thus, one gets the standard pairing Hamiltonian: 

H P = Z^Pj-sL^fi A l = L «%flj-m ( 6 ) 

j jf m>0 



This 2-body Hamiltonian, however, is exactly solvable even for systems with more 
than one pair since it can be viewed as Richardson-Gaudin model |24|. For example, 
the relevant equations for the proton-neutron T — 1 pairing that were given in Ref. [23) 
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as well as by Links et al. (25], and Asorey et al. (26] are: 



L r\ Mo M — 1 1 

m2et-v a p%v a -v p j^wy-va 



M [ M-T -y M 



*=E 



— V a - Wy 5 — Wy - W S —j " ' 

The spectral parameters v a have the same meaning as pair energies. The Wy parameters 
are related to the iso-spin symmetry of the proton-neutron pairing. By drooping the 
terms that contain the Wy parameters one arrives at the Richardson's exactly solvable 
pairing for one type of particles (27]. If one considers only one pair case (M — 1 = p) 
in L single particle levels with double degeneracy of each single particle level i (£2, = 
(2j'+l)/2 = l) then one has the one pair energy eigenvalues of the EPI |5] and the 
standard pairing (6]: 



E=z, 



1 £ 1 



-\2Ei-Z 



This is a special case of the (p = 1) solution for the extended pairing model (9): 

E\ = zf -G(p- 1), i = £ — = £ 2e,„. (8) 

° ii...i p e 'ii-ip z n=l 

In the above equations, we intentionally kept the notation for the 2-body pairing and 
the A-body pairing slightly different to emphasize their different structures. 



5.1 Binding Energy of the Sn and Pb Isotope Chains 

Deformation is common in very heavy nuclei and this often justifies the success and 
application of the Nilsson model. For the purpose of our model, we use deformation 
parameters from Ref. [28] and experimental binding energies from Ref. [29|. Theo- 
retical relative binding energies (RBE) are calculated relative to a specific core, 152 Yb, 
100 Sn, and 208 Pb for the cases considered. The RBE of the nucleus next to the core 
is used to determine an energy scale for the Nilsson single -particle energies. For an 
even number of neutrons, we considered only pairs of particles. For an odd number of 
neutrons, we apply Pauli blocking to the Fermi level of the last unpaired fermion and 
consider the remaining fermions as if they are an even fermion system. The valence 
model-space consists of the neutron single-particle levels between two closed shells 
with magic numbers 50-82 and 82-126. By using d§), values of G are determined so 
that the experimental and theoretical RBE match exactly. 

Here we discuss mostly the Sn isotopes since the Pb and Yb isotopes were discussed 
in more details in Ref. (TT[ and Ref. (TO] , In Figure [5] are shown the results for Sn as 
calculated by using the lt)u Sn as core and zero RBE nucleus. The single-particle energy 
scale is set by the binding energy of 101 Sn. The inset shows the fit to values of G that 
reproduces the experimental data exactly. The two fitting functions are: ln(G(A)) = 
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365.0584 -6.4836A + 0.0284A 2 and ln(G(A)) = 398.2277 -7.0349A + 0.0307A 2 for 
even/odd values of A. The solid line gives the theoretical RBE of the Sn isotopes using 
these fitting functions. 




The Sn isotope chain is unique in the sense that there are two doubly magic mem- 
bers, the Sn and Sn, that allows us to use Sn as zero RBE system as well. 
In Figure [6] are shown the results for Sn isotopes when using 132 Sn as zero RBE sys- 
tem. In this case, there are again good even/odd quadratic dependence of the ln(G(A)), 
however, as for Pb case ]TT) there is a simpler expression that works for even and odd 
systems simultaneously. In this case we have G(A) = adim(A)-^ with a = 259.436 
and = 0.9985. 



6 Conclusion 

In this paper we have presented evidence for the need to use accurate proton and neu- 
tron masses, or at least a properly weighted nucleon mass in order to improve on 
the accuracy of the binding energy of light nuclei as computed with the next generation 
computer codes. This will also allow better understanding of the NNN-, NNNN-, and 
A-body interactions in nuclei either derived from ChPT or from a phenomenological 
considerations. Therefore, one has to build A-body computational technology in the 
next generations of nuclear modeling codes. 

While the motivation for considering A-body interaction in the light-nuclei is strong 
as based on the ChPT QCD derived interactions, one is left to wonder if A-body in- 
teractions are also relevant to heavy nuclei. The results obtained with the help of the 
Extended Pairing Interaction, in particular the Sn isotopes discussed here, seem to 
confirm the idea that A-body interactions are needed to understand better the binding 
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energy of heavy nuclei. Often the imagination cannot capture all the possible implica- 
tions and uses of an exactly solvable model. Beside the current applications of the EPI, 
one can also see that it would be a useful verification tool for A-body computational 
codes as well. 
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